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Abstract 

Instanton-like transitions ("shadow processes") are considered in (1+1) di- 
mensional models with one scalar field whose potential is a quadratic well with 
a cliff. The corresponding classical boundary value problem is solved, and the 
semiclassical transition probabilities are found in a rather wide range of energies 
and the numbers of initial particles. 
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1 Introduction 



Much effort is being made to understand the high energy behavior of the amphtudes 
of instanton-like processes. Being suppressed at zero energy by the factor e~"^", where 
5*0 is the action of the instanton, these amphtudes exhibit the exponential growth with 
energy, and, in the leading order of the perturbation theory around the instanton, 
cross the unitarity bound, typically at the sphaleron mass scale . This observation 
has lead to an exciting speculation that the electroweak baryon-number violating 
processes can become observable in the TeV energy region H (for reviews see 
refs.ii). 

There are strong arguments implying that the behavior of the instanton-like cross 
section is exponential. 



a = exp 



9' Uo 



where the energy scale Eo is of the order of the sphaleron mass in the electroweak 
theory and many other models; the function F is the central object of current studies. 
Until now, only low energy asymptotics of F{E/Eq) is understood, and a few terms 
of the corresponding series have been explicitly calculated within the perturbation 
theory around the instanton ||, |3, ||, If. The pecuhar form of the cross section 
suggests that even at energies comparable to the sphaleron mass, where the function 
F cannot be found by the perturbation theory around the instanton, there might 
exist a semiclassical-type procedure for calculating the amplitudes. The key problem 
for such a procedure is the non-semiclassical nature of the two initial hard particles 

Recently, several approaches that may allow to overcome this problem have been 



proposed [|Tl|, |T2|, |T3|, |T^, 0. One of them [|TT], |T2| is as follows. Instead of considering 
two-particle scattering, one examines the maximum probability of transitions among 
all initial states with given energy E and number of initial particles n = u/g'^, where 
u is fixed in the limit g"^ —>■ 0. In this regime, the total probability indeed has the 
exponential form (|l|), where the function F now depends on u. At finite u, there exists 
an entirely semiclassical procedure for calculating the function F{E / Eq,v) which, 
in its final version, is reduced to a well-defined classical boundary value problem 
[|l^. Though not a priori justified, a rather natural assumption is that the limit 
z/ — s> of the function F{E/Eq, v) is smooth and its limiting value F{E/Eq, z/ — * 0) 
gives at least the upper bound for, and may even coincide with the two-particle 
function F{E/Eq). The latter conjecture has found substantial support from high- 
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order perturbative calculations around the instanton [|T7| . 

The prescription for the evaluation of the function F{E / Eq^v) is as follows (see 
ref.pBI for details). One searches for a solution to the field equations, 



on the contour ABCD in the complex time plane which is shown in fig.l, with the 
following boundary conditions, 

ct>{k) = ^(/ke-^-'=*' + e'f^e^-"'') when t' ^ -oo, 

0(fc) = ^^(6ke-*^'=* + 6*„ke'""=*) when t +00, (3) 

where T and 6 are two free parameters of the field configuration which are related 
to the energy and the number of initial particles, t' = t — iT/2 is the real parameter 
on the part AB of the contour, and /k and 6k are arbitrary functions of the spatial 
momentum k. In other words, the boundary conditions @ require that the ratio 
between the negative- and positive-frequency parts of the asymptotics of the field 
is equal to e~^ in the initial state, and 1 in the final state. The solution to this 
boundary value problem is complex on the contour of fig.l (except for = 0), but is 
real on the part CD. It describes the most probable "microcanonical" process among 
all initial states with energy 

E = J dkujke'^f^fk = J dkukb^bk 

and number of particles 



n = J dke'flf^. 



The classical action of the field evaluated along the contour of fig.l, which we will 
denote by iS (so that 5* is real and positive for purely Euclidean fields, such as an 
instanton), satisfies the following relations, 

dS ^ dS 

dT ' 89 ' ^ ' 

and the function F is merely the Legendre transform of the double action, 

^F(E,n) = 2S - ET -n9. (5) 
9 
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Obviously, it satisfies the following relations, 

It is likely that the above boundary value problem, because of its complexity, can 
be solved only by numerical methods, except for the region of low energies, E <^ Eq, 
where F can be found by the perturbation theory around the instanton. However, for 
the qualitative understanding of the structure of the solution as well as the behavior 
of the function F{E/Eo, z/), it is desirable to invent a simple model where the solution 
and the function F can be found analytically, at least in some non-trivial range of E 
and n, by techniques different from the usual perturbative expansion around a single 
instanton. 

In this paper we present such a model. We consider a (l+l)-dimensional theory of 
one scalar field, whose potential has a local minimum at = 0, but is unbounded from 
below. The decay of the false vacuum = is described by a bounce configuration 



IS], which shares many properties of instantons in gauge theories. Making use of 



this bounce solution, one can consider a "shadow process" [|19|, 0, i.e. scattering of 
particles through the formation of a bubble. The shadow process is a direct analog 
of instanton-induced scattering and it can be treated by the same methods as those 
used in gauge theories. Though not precisely defined at high energies, the "shadow 
process" is a good choice for testing the semiclassical technique proposed in ref . ||16[| . 

In the model that we will study in this paper, there exists an additional large pa- 
rameter besides the inverse coupling constant. When this parameter tends to infinity, 
the potential has the form of a quadratic well with a cliff. In this limit, the boundary 
value problem is solvable in a rather wide range of E/Egph and u, and the solution 
is a non-trivial generalization of the dilute instanton gas. The specific models to be 
considered in this paper are described by the exponential interaction, exp(A0), and 
power-like one, 0^, where A and are large. 

For the two models, the results are qualitatively the same. The periodic instanton, 
which is the solution to the boundary value problem at 6 = 0, can be found in the 
whole range of energies where it exists, i.e., < E < Egph- For u of order 1 the 
function F can be found for all energies smaller than some critical value Ecntiv) 
where the exponential suppression is reduced by a factor of order A^^ ^ 1 (i.e., at 
the critical energy F ~ Sq/X)\ the smaller the number of initial particles, the larger 
the critical energy. For even smaller z/, namely, v ~ A^^, the approximations made in 
this paper fail at relatively low energies, where the exponential suppression is reduced 
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by a factor of order one only. The limit z/ ^ is smooth and in the exponential model 
the function F{E/Eq) coincides with the naively extrapolated one-instanton result, 
that shows that all higher terms of the perturbative expansion around the instanton 
vanish. In the power model, the result does not coincide with the one-instanton 
formula and corresponds to summing up the whole perturbation series, where each 
term is calculated to the leading order in 

The fact that the function F{E/Eq,u) can be calculated beyond the perturba- 
tion theory is a remarkable feature of our model. Unfortunately, our results for the 
most interesting case u ^ are rather limited (we can only trust the reduction of 
the exponent by a factor 2/3 and e"^/^ in the exponential model and power model, 
respectively), so they do not tell much on whether the two-particle instanton-like 
scattering is strong or not at high energies. 

One question that can be addressed quantitatively in these models is whether 
multi-instanton processes become important at those energies when one-instanton 
amplitudes are still exponentially small. This conjecture pT| , that has lead to 



a hypothesis of premature unitarization, is discussed in this paper within the expo- 
nential model. For processes with relatively large number of initial particles, which 
we can safely treat until the exponential suppression is strongly reduced (i.e. until 
E = Ecritiy)), we find that multi-instanton contributions are more suppressed than 
single-instanton ones at all energies up to Ectu- We think this is a strong argument 
against the premature unitarization hypothesis. 

The paper is organized as follows. Sects. 2-5 are devoted to the exponential 
model. In Sect. 2 we describe the exponential model and find the analogs of the 
sphaleron and the instanton. Sect. 3 is the central part of the paper and contains the 
construction of what we call the "improved dilute instanton gas approximation" . The 
solution to the boundary value problem is found, and the function F is calculated 
for the two cases z/ ~ 1 and v ~ 1/A. The periodic instanton solution and multi- 
instanton contributions are also discussed. In Sect. 4 we try to go beyond the improved 
dilute instanton gas and describe the periodic instanton up to Egph] unfortunately, we 
are unable to obtain essentially new information on the processes with small number 
of initial particles at high energies beyond the improved dilute gas. Sect. 5 is devoted 
to a related problem, namely, induced vacuum decay, which can be also considered 
as a classical boundary value problem. The exponential suppression of the vacuum 
decay is not reduced when the number of initial particles is small, up to the energies 
where our approximation breaks down. The model with 0^ interaction with large 
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A^, where our approach requires some modifications, is considered in Sect. 6. Sect. 7 
contains concluding remarks. 



2 The model 



2.1 Shadow process 

The model to be considered in this paper, except for Sect. 6, is a (1 + 1) dimensional 
theory containing one scalar field 0, with the lagrangian 

1 



L 



.(d,4>y-v{ct>). 



The potential V"(0) has a mass term and an exponential interaction term. 



v{4>) 



m 



2 1 

m V 



exp 




(7) 



Note that the interaction term has negative sign, so the potential is unbounded from 
below at large positive 0. The theory contains two massless parameters, v and A. The 
parameter 1/f is the coupling constant that governs the perturbation theory. The 
probabilities of semiclassical processes we are interested in are exponential in v. in 
general, expressions like eq.(|l|) will have the form 



exp 



-F V 



2 ^ X 



exp 



(i.e., Ino" oc —v"^ at fixed E/Eq and z/), where Eq oc and the number of initial 
particles is n = v'^u. In this paper we concentrate on the exponent F and make no 
attempts to estimate the pre-exponential factors. For the semiclassical expansion to 
be consistent, we take A -C f . On the other hand, A itself needs not be small, and we 
assume throughout this paper that A ^ 1. So, we take 

w > A > 1. 

The behavior of the potential is shown in fig.2. Due to the large value of the 
parameter A, the potential has the form of a quadratic well with a cliff. At (p = v, the 
mass term and the interaction term in the potential cancel each other, so V{v) = 0. 
When (f) decreases, the exponential term rapidly decays and becomes negligible as 
compared to the mass term, so for (t> — 0) ^ v/\ the potential reduces to its quadratic 
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part, m^0^/2. When is larger than the interaction term, however, becomes 
dominant, and the potential falls down rapidly. The potential has a local minimum, 
which at large A is placed almost exactly at = 0. This minimum corresponds to the 
false vacuum, whereas the true vacuum is = +oo. 

The decay of the false vacuum is a non-perturbative phenomenon with the ex- 
ponentially small amplitude, which is typical for quantum tunneling processes. This 
decay can be described semiclassically by the bounce configuration, which is a solution 
to the Euclidean field equation with finite action, that obeys the boundary condition 
0(oo) = 0. The action of the bounce determines the exponential suppression factor 
for the false vacuum decay. 



Since the bounce is the analog of the instanton, it is suggested in refs.[jT9|, ^ that 



theories with unstable vacua may be considered as models for studying instanton- 
like scattering. Instead of transitions from one vacuum to another induced by particle 
collisions, one studies "shadow processes", where both initial and final states refer 
to the false vacuum = but intermediate states contain a bubble of large positive 
field (0 > f in our case). The simplest shadow process is one described by the bounce 
itself, which corresponds to the transition at zero energy from the false vacuum (0 = 
at r = +oo) through some bubble-type configuration (0(x) at r = 0) back to the 
false vacuum (0(a;, — oo) = 0). We note in passing that the very definition of the 
shadow processes is of semiclassical nature, and the notion of the shadow processes 
may become ambiguous in high energy scattering. 

Along with the bounce, there exists another configuration relevant to shadow 
processes, namely, the critical bubble, which is a static and unstable solution to the 
field equation with finite energy. The physical significance of the critical bubble is 
that it is the minimum static energy configuration on the border, in the configuration 
space, between the false and true vacua: being slightly perturbed, the critical bubble 
either evolves (in real time) into an expanding domain of a new phase that eventually 
fills up the whole space, or shrinks down to = 0. The critical bubble is the analog 
of the sphaleron in gauge theories (in particular, its free energy determines the rate 
of the false vacuum decay at high enough temperatures). In many theories (but not 
in our model, see below), the sphaleron mass determines the energy scale where the 
instanton-like amplitudes naively become large. For convenience, we will use the 
gauge theory terminology and call the critical bubble and bounce by the sphaleron 
and instanton, respectively. 
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2.2 Sphaleron 

The sphaleron is the static solution to the field equation, 



d^cf) = w^cf) — Xrn^v exp 



This equation has the same form as the equation of motion for a classical particle 
in the upside-down potential —V{(f)) with x playing the role of time. The sphaleron 
corresponds to the motion of the particle from = Oatx = — ooto0 = fata; = O 
where it recoils from the wall and goes back to = at x = oo. At large A, when the 
wall is steep, the recoil occurs in a small interval of "time", and most of the "time" 
the particle moves in the region where the nonlinear term in eq.(P) can be neglected 
and the potential is quadratic. So, outside the nonlinearity region, the sphaleron 
configuration is 

= e-'"l"l. (9) 

One can check that the region of nonlinearity, where the exponential term in eq.(^ 
is comparable with the other terms, is of order 

1 

Am ' 

which is much smaller than m~^. The energy of the sphaleron, 

E = Jdx [^(9,0)' + ^0) 

is then dominated by the contribution from the region where the field is linear, x ^ 
(Am)~^; the contribution of the region of nonlinearity is suppressed by a factor 1/A. 
A straightforward calculation gives 

Esph = mv^ (10) 

in the leading order in 1/A. 

The fact that the sphaleron mass comes mostly from the linear region is a specific 
feature of our model, which allows, in particular, to calculate the momentum distri- 
bution and the total number of particles emitted by a sphaleron decaying back into 
the false vacuum. The field of the decaying sphaleron is linear, 

/d Jc 
^ ^ (&fce~^"^-*+^^" + 6*_^e''"'=*-''="), 
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and, if the sphaleron begins to decay at t = 0, obeys the initial conditions t 
0) = 4>sph{x), (j){x,t = 0) = 0. One finds from eq.(||) 

1 mv 



The total number of particles emitted by the decaying sphaleron is 

2 

nsph = I dkblbk = -v^ ■ fH' 

TT 



As expected, the number of particles is of order v 



2 



2.3 Inst ant on 

The instanton is an 0(2) symmetric solution to the Euclidean field equation, 

o2 , , 1 <90 2 , X 2 

a.^(p = + -T^ = m <p- Am v exp 




(12) 



= (r, x), r = Vt^ + X 



The appropriate mechanical analog is now a particle moving in the potential —V{(f)) 
with the friction coefficient inversely proportional to "time" r. The motion of the 
particle begins at some 0(0) > v with zero velocity and ends at = at infinitely 
large "time". The initial potential energy is lost at r = oo because of friction. 

We are unable to solve eq. ([1^) analytically. Nevertheless, by making the following 
conjecture on the structure of the field configuration we are able to obtain the solution 
to the leading order in 1/A. Namely, in analogy with the sphaleron, we assume that 
the size of the region of the Euclidean space-time where the field is nonlinear, ro, is 
much smaller than the inverse mass, ro <^ m^^ . Outside this nonlinear "core" the 
field is free, i.e. the nonlinear term in eq.(|l^) is negligible as compared to the kinetic 
and mass terms. This outer region r ^ ro can be in its turn divided into two regions: 
ro <^ r ^ mT^ where the mass is unimportant and the field is both free and massless, 
and r > m'^ where the field obeys the ordinary massive Klein-Gordon equation. 
Shortly speaking, we assume that the field shows different behavior in the following 
three regions. 



(1) r < ro, where d'lcj) = — Am^f exp 



2A ^ 



(2) ro ^ r ^ m ^, where d'lcj) = 0, and 
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(3) r ^ m ^, where 9^0 = m^0. 

In the region 2 we will match the two solutions found in the regions 1 and 3 thus ob- 
taining the complete field configuration. We will justify our assumptions a posteriori, 
and now we construct the solution explicitly. 

First, consider the regions 1 and 2, where r <^ and the mass term can be 
neglected. Since there this no mass scale left over, the equation is scale invariant. In 
fact, it is the Liouville equation, which in terms of the new variable 

X = --0 (13) 

V 

has the form 

dlx = (14) 

where 

2e^ . X 

At the moment we have to consider the 0(2) symmetric solution to the Liouville 
equation. 



X(r)=ln c+— , (16) 

y Cd j 

where c is an arbitrary constant. Thus in the regions 1 and 2 the field configuration 
is 

where c is yet to be determined. Notice that at r ^ ca, where 0(r) ~ Inr, the field 
is free and massless. This means that the radius tq separating the regions 1 and 2 is 

To = ca. (18) 

In the regions 2 and 3, i.e. at r ^ Tq, the field is linear, (9^ — m^)0 = 0, so the 
configuration is given by the 0(2) symmetric solution to the Klein-Gordon equation, 

0(r) = a_K'o(mr), (19) 

where a is some constant and Kq is the modified Bessel function. If rg ^ m^^, we 
can match the two solutions, eqs. ([17|) and (|19D, in the region 2. We write 

A V cd^ , 





mr 


= —a 


In h 7 

2 . 
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where we made use that at small (mr), Ko{mr) = — ln(mr/2) —7, where 7 = 0.577 . . . 
is the Euler constant. Taking into account eq. (|T5|) , we find the parameters a and c, 

c = AV2("+^). (20) 



So, in the regions 1 and 2, the solution is given by eq.(|T^, and in the regions 2 and 
3 it is determined by eq.(|l9D, with the values of a and c fixed by eq . (|20|) . 



It is now straightforward to justify our basic assumption that Tq -C m ^. We find 
from eq.(|l8D 

ro ~ Ae~^— (21) 
m 

which is smaller than by an exponential factor e~'^. So our basic assumptions 
about the structure of the field are indeed valid with exponential accuracy. 

Thus, we have found the instanton by solving separately the field equation in the 
regions of small and large r and matching the two solutions at intermediate r. The 
assumptions made on the structure of the instanton solution are self-consistent at 
large A. The field configuration consists of an exponentially small core, where is 
nonlinear, and a large tail that extends from rg ~ e~^m~^ to . In fact, it is 
straightforward to see that in the whole Euclidean space-time our solution satisfies 
the field equation with the accuracy e~'^. 

Note that the value of the field at the center of the instanton is finite in the 
limit A ^ 00, 

0(0) = -^lnc = 2t;. (22) 
A 

The instanton action is dominated by the contribution of the kinetic term in the 
lagrangian, and comes mostly from the region 2, 

5o = 2vr I rdr\{drct^f = ^ (l - ^ + 0{\-')^ . (23) 

Recall that we take f ^ A ^ 1, so the vacuum decay probability is exponentially 
small. 

2.4 Leading order instanton approximation 

In the leading order of the perturbation theory around the instanton, the two-particle 
cross section at energy E can be found by the technique of ref. One writes for the 
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exponent of eq.(|l|), 

F{E) = 2So-ET- J dkR*{k)R{k)e-'^'=^, (24) 
where R{k) is the Fourier component of the asymptotics of the instanton field, 

0(0;, t)= [ J^l_ i?(fc)e-^-+'^-, r oo, (25) 

and T is related to the energy by the following condition, 

E = J dkuJkR*{k)R{k)e~'^'^^. (26) 



The Fourier components in eq.(p^ can be calculated from eqs.([T9|) and (^ 



The integrals in eqs.(|^) and (p6D are then straightforward to evaluate, and for E ^ 
Egph one obtains the leading order instanton approximation result. 

Naively, eq.(P7D implies that the cross section becomes large at 

E^Eo = e^E^pf,. (28) 
This scale is much higher than Esph, which is a peculiar feature of our model. 

3 Improved dilute instanton gas approximation 

3.1 Difficulties of the ordinary dilute instanton gas 

At low energies, the classical boundary value problem, eqs.(|^) and (0), can be solved 
by making use of the dilute instanton gas approximation ^ . Typically, the dilute 
instanton gas approximation breaks down at ~ Esph', this is the case, for example, 
in the (l+l)-dimensional Abelian Higgs model []TB[ and the electroweak theory 



Let us see, however, that our model is peculiar in this respect: the ordinary dilute 
instanton gas approximation breaks down at energies much smaller than Egph- 

Let us consider, for example, the boundary value problem at 6' = 0. In that 
case the solution is periodic in Euclidean time with the period T |2^. The basic 
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assumption of the dilute instanton gas approximation is that the Euchdean solution 
may be represented by the sum of widely separated instantons, which, to the leading 
order, do not distort each other's interiors. Thus, the periodic instanton field is 
approximated by the sum 



I 



'inst 



X 



Xl), 



(29) 



where, by periodicity, xi = {IT, 0) are the positions of the instantons in the chain 
(there are no anti-instantons in our model). 

At first sight, the dilute instanton gas would be a good approximation if T is much 
larger than the core size given by eq. (|2T|) . However, we will see immediately that 
this condition is not sufficient for a two-fold reason. First, the field in the center of the 
instanton, eg. (p^) , is strong enough, so that the exponential interaction is operative. 
Therefore, even relatively small distortion of the field due to other instantons may 
have strong effect on the instanton interior. Second, the field of each instanton falls 
off rather slowly at |x — x/| ^ m~^, so that neighboring instantons may produce a 
collective effect on a given instanton. 

To establish the region of validity of the dilute gas approximation, let us consider 
the interior region of an instanton sitting at x = and write explicitly its contribution 
to the sum in eq. (|29|) , 

0(x) = (pinst{x) + 0(x), (30) 



where 



4){x) 



Hnst 



[x - Xi) 



is the distortion of the field due to other instantons. For the dilute instanton gas to 
work, the distortion (p should not change the field equation considerably. In particular, 
at a; = one should have 



exp 



exp 




- 1 



which means 



Recalling that, due to eqs.(|l9D and 



(x = 0) < 



2v 



(31) 



ct>{x = Q)=Y.-K^{m\l\T), 

1^0 ^ 
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one finds that eq.(|3T|) is valid only when 

T > m-\ 

Obviously, this condition is much stronger than the naive estimate T ^ tq. We will 
see in Sect 3.2 that the corresponding condition on energy is 



so that the ordinary dilute instanton gas approximation breaks down at energies much 
lower than the sphaleron mass. 



3.2 Improved dilute gas 

The ordinary dilute instanton gas approximation breaks down long before the in- 
stanton cores begin to overlap. One may expect, however, that the picture of well 
separated instanton-like objects is valid in a wider region of energies, but each instan- 
ton is strongly affected by the collective field of the others. When the instanton core 
is much smaller than the separation between the instantons, this collective field is 
constant in the instanton interior, so the distorted instanton is still 0(2) symmetric. 
Furthermore, outside the instanton cores, the field is still linear, so that we are still 
able to make use of the assumptions described in Sect. 2. 3. 

To obtain the solution to the boundary value problem, eqs.(0) and (|), we begin 
with the region outside the cores. Let the distorted instantons be located at a; = and 
T = (±Ti ± IT, 0), / = 0, 1, 2 . . ., where Ti is yet unknown parameter (see fig. 3). At 
6^0, only two of these instantons (with / = 0) have unit strength, while others are 
fake instantons with strength exp(— (cf. ref.[|16|). Recall that we are interested 
in the field on the contour ACD in the complex time plane shown in fig. 3, so we are 
able to ignore interiors of all instantons except for one sitting at r = Ti. 

Thus, the field outside the core is 

0(x, r) = a 2" e^''''^ Ko (m^{T - Ti - IT f + + (m^J{T + Ti-lTy + x^ 

l=-oo ^ ^ ^ ^ 

(32) 

where a is yet unknown constant that may be different from eq.(|20D. The field (|32D 
automatically satisfies the boundary conditions, eq.(|^). The Fourier components of 
the asymptotics of the field in the initial and final states are 

— 2e— cosh(cjfcTi) 



fk - fk 




1 - e-'^kT- 
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The energy E and the number of initial particles n are related to the parameters T, 
Ti and 6 as follows, 



E = y dkuke'flh = Tia' J dk——^ 



2 cosh cjfcTi 



2 / (e"*^*"^^ + e~'^*('^~'^^^^^ 



dkukblh = TTa J dk _ ^^^^T-ey ' (^4) 

n = jdke fj, = na J - _ • (35) 

In fact, these three relations are consequences of eq.(^ and the condition of equilib- 
rium of forces acting on the instanton sitting at r = Ti. Eqs.(0) and (|35|) enable one 
to express T, Ti and 6 through the energy and number of initial particles, once the 
value of a is known. Alternatively, one may consider the boundary value problem at 
given values of T and 9, and relate these parameters to E and n at the very last step. 



In that case eq.(^J) should be used to express Ti in terms of T and 6. 

To find the configuration inside the core, let us first consider the region close to 
the instanton located at (Ti, 0), but still outside the core. In this region, all terms in 
the sum (^) are constant with the accuracy tq/T, except for one that refers to the 
instanton (Ti, 0). Eq. (p^) then reduces to a form that is 0(2) symmetric with respect 
to the point (Ti, 0), 



(x, t) = a 



(36) 



Ko(Tn^ir-T,y + x^j+f^, 
where / is a constant that depends on T and 9, 

oo oo 

/ = 2^e-'^ir(nmT) + ^ e-l'l^i^o(m|ZT - 2ri|). (37) 

1=1 l=—oo 

A convenient representation for this constant is 

_ f dk 2e-"^^-^ + + 
^"7 2^ 1 - e--^-^ • 

Eq.(^) can be considered as the boundary condition for the Liouville equation 
that determines the field inside the core. Since this boundary condition is 0(2) 
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symmetric, the field inside the core is still given by the 0(2) symmetric solution of 
the Liouville equation (p!4D, 

^''^ + ^-1, (39) 




9 I ' 

ca^ J 

where a is given by eq. (|15]) and c is some constant that can differ from that of the single 
instanton, eq.fpOl). If the core size, which is Tq = ca, is much smaller than both Ti and 
T/2 — Ti, one can match eqs.(^) and (^) in the region Tq ^ (r — Ti)^ + x"^ <^ Tf , 
(T/2 — Ti)^ where is free and massless. In this way we obtain the values of the 
parameters c and a, 

c = A^e-^^+i-fl, (40) 

(41) 

Eqs.(|36|) - ( PD determine the solution to the boundary value problem. For given 
T and 6, the value of Ti is determined by eg . (^If) . The value of / is then given by 



eq . (|38|) , so that all parameters become known. Outside the cores, the field is exactly 
the same as in the ordinary dilute instanton gas approximation: the configuration 
( |32D with a given by eq.pl|) is precisely the sum of undistorted instanton fields. On 
the other hand, the size of the core and the field inside the core are modified. Namely, 
the core size is 

ro = ca = 2Aexp(-A + 2/-27) — , (42) 

m 

while the field configuration inside the core is the 0(2) symmetric solution of the 
Liouville equation with the parameter c given by eq.(^). 

We will call this procedure the "improved dilute instanton gas approximation". 
This approximation is valid provided the core size, eq. (^2]) , is much smaller than Ti 
and T/2 - Ti, 

Aexp(-A + 2/)l<T, (T/2-Ti). 
m 

The explicit formulas will be given in Sects. 3.3 and 3.4 where we will also discuss the 
actual region of validity of the improved dilute instanton gas approximation. We note 
in passing that at T ^ m~^, the improved and ordinary dilute gas approximations 
coincide. Indeed, in that case all terms on the right hand side of eq.(|37D are much 
smaller than 1, so that / ^ 1. The value of the parameter c, eq . (|40|) , is then the 
same as for an isolated instanton, eq.(^Ol), and the field in the instanton interior is 
precisely the field of an isolated instanton. 
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When 6 is fixed, and T decreases, the value of / increases. Since / enters the 
expression for tq, eg . (|4^) , in the exponent, the core size increases rapidly when T 
falls below m"^. At some T = T^ru which depends on 6, the core size becomes 
comparable to the distance between the instantons and our approximation breaks 
down. We will see that this happens at energies comparable to, or even much larger 
than the sphaleron mass, depending on the value of 6. 

Now let us evaluate the action for our configuration. Since the core is exponentially 
small, the action comes entirely from the region where the field is linear, so one has 
(recall that we choose the convention that the action is real and positive for real 
Euclidean fields) 

S = --j dtdx{{d,(t>f - m^^), (43) 

where the integration over time is performed along the contour of fig. 3. This contour 
can be divided into two parts: {iT/2 — oo,iTi), and (iTi, +oo). In the first part the 
field is given by the initial state asymptotics, while in the second part it is given by 
the final state asymptotics. Substituting these asymptotics, eq.(H), into the integral 
(^31) , one obtains 

S=\j dk{fir_,e-^^^-'^^^^'<^ - hUe-^^^-'^^^) 

The explicit form of the Fourier components is given by eg. (|33|) , so we find 

_ 47rf 2 r dk 1+ e'^kT-e _^ ^-ukT-e-\-2wkTi _^ g-2<^fcTi 
J 2^k 1 _ e— fe^-e • 

Comparing this expression with eg . (^8|) , we obtain 

The integral in this eguation diverges logarithmically in the ultraviolet region. How- 



ever, we have to cut this integral off at the scale I/tq, since eg . (|33|) is valid only at 
momenta much lower than I/tq. So, we have 



Recalling that the core size tq is related to / by eg. (021), we obtain finally 



S^'-^(,-L^0('^)] (44) 



A I A I A 
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It is straightforward to verify that this action as function of T and 9 satisfies eq.(^. 

To obtain the exphcit formulas for the transition probabihties at a given number 
of initial particles, we have to consider various limiting cases. Since there exists an 
extra large parameter A besides the "coupling constant" f , we distinguish the cases 
V = n/v"^ 1 and z/ ~ 1/A where n is the number of the initial particles. We will see 
in what follows that these cases correspond to ~ 1/A and 6' ~ 1, respectively. Let 
us discuss the two regimes in turn. 



3.3 

Let us first consider the case of small 6, 6 ^ 1/X. At small 6, the condition of energy 
conservation, eg. (P^) , implies that Ti ^ T/4. Let us also assume that T ~ (Am)~^; 
we will see that this regime corresponds to ~ Esph- For these values of T and 6, 
the number of relevant terms in eq. f^Tl) is of order A, which is large. In other words, 
there are many fake instantons that determine the solution on the contour of fig.3. 
For the actual calculation of /, it is more convenient to make use of the integral 
representation, eq. (|38D . At A; <^ Am, the integrand in eq.(^8|) can be expanded in the 
following way. 



So, one has 



^ r dk 4 r dk 



- 1 + 0(A- 



2ujk uJkT + 9 J 2ujk 
The second integral on the right side diverges logarithmically in the ultraviolet, and 
should be cut-off at A; ~ Am. We obtain 

f = ^h{9/mT)-\n\ + 0{l), (45) 

where 

, , 2 arccosh/3 

JiKP) = ,n-2 1 • 

71 ^/I3^ - 1 

Note that by our assumptions, 9 ~ 1/A, mT ~ 1/A, the first term in eq.(|45|) is of 
order A. Analogously, the energy and the number of initial particles are obtained 
from eqs.(P^) and (P3|), 

E = (^) ' /2 {9/mT) E,,,, n = ' f, (O/mT) n,,,, (46) 
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where 

arccosh/3 




(/32 - 1)3/2 J ' 

/3(/3) = ^(/l(/3)-/2(/3)), 

and the energy and the number of particles for the sphaleron are given by eqs.([TO|) 
and ([Tl|). We see that the regime 6* ~ 1/A, mT ~ 1/A occurs when the energy E and 
the number of particles, n, are of the same order of magnitude as the corresponding 
quantities for the sphaleron, Egph and n^pft. Since risph ~ v^, we conclude that the 
regime ^ ~ 1/A corresponds to u = n/v"^ ^ 1. 

The exponent for the total probability is then straightforward to evaluate, 

_ Snv^ / Atx 

^ ~ir \ " A^" 

The condition for the validity of the improved dilute instanton gas approximation is 
ro -C T. Since T ~ (Am)~\ and 

1 



l-——M9/mT)). (47) 



the approximation is valid if 



ro~Aexp(-A + 2/)- 
m 



;A (e/mT) » \. (4J 



XmT A 

Comparing eqs.(^) and (^Sp we see that we can trust the dilute gas approximation 
up until the suppression becomes much weaker than the instanton suppression, 

A 

In the particular case ^ = 0, the field configuration is the periodic instanton that 
describes the most probable process at given energy E [2^ (the maximization over 
all possible values of the number of initial particles n can be seen from the relation 

= = 0). Since the functions /i, /2, /a are normalized in such a way that 

/i(0) = /2(0) = /3(0) = 1, we have for the periodic instanton 

E _ / 47r \2 
Esph \\mT ) 

and the function F{E) can be found explicitly. 
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The exponent for the probabihty, {—F{E)) for the periodic instanton is shown in 
fig.4. It increases from — 25*0 at low energies to zero at E = Egph- The number of 
particles for the periodic instanton increases with energy, 

E 

^sph 

SO that the average energy per particle remains constant. The core size of an individual 
instanton rapidly increases with energy, 



'"o ~ -; — exp 
Am 



-All ■ ^ 



E 



sph ^ 



Strictly speaking, there is an energy region close to the sphaleron mass where the 
improved dilute instanton gas approximation does not work. According to eq.(^), 
this region is 

1 ^ ^ 
i ~ — . 

Esph A 

We will discuss this region in detail in Sect. 4 and show that the behavior (|49|) indeed 
persists up to Egp^. 

Let us now consider the general case ^ 7^ 0, still assuming ^ ~ 1/A. To find the 
probability of the transition at fixed energy E and number of initial particles ra, which 
are of the order of Egph and Usph, respectively, one finds T and from eq. (|46|) and then 
substitutes them into eq.(^). The result is presented in fig.5. For a given n < nsph 
the function F{E) starts from E = 7rnm/2, where ^ = and the configuration is just 
the periodic instanton, and reaches zero at some Ecru determined by the condition 

^ /i (O/mT) = 0. (50) 



XmT 

Very close to Ecru, however, the improved dilute instanton gas approximation is not 
reliable, see eq.(^). 

In the whole energy interval, F is a monotonically increasing function of energy. 
As in the case of the periodic instanton, at the point where the approximation breaks 
down, the function F is of order Sq/X. The critical energy at which F becomes small, 
is a function of the number of initial particles: the smaller the number of initial 
particles, the larger the critical energy. 

At 1/A <^ n/usph <^ 1, the critical energy can be calculated analytically. In that 
case one has 6/mT ^ 1, so that one makes use of the asymptotics of /i, /2 and at 
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large values of their argument, 

21n/5 



fm = 



vr (3 ' 
_ 2 1 

ln/3 



It is then straightforward to express T and 9 through E and n from eg . (^Bf ) and then 
obtain the critical energy at which eq. (|50D is satisfied and the instanton suppression 
is strongly reduced. One finds the following exponential dependence, 

Ecrit{n) (T^'^nsph\ 
- exp ' ' 



Esph \A n J 

Note that we consider n/usph ^ 1/A, so Ecnt is still much smaller than the energy 
scale where the two-particle amplitudes naively become large, eq.(^). 



3.4 ^ - 1 

In the case 6* ~ 1, the sum in eq.(P7D is saturated by a finite number of terms. Let 
us consider the case mT -C 1, which is, of course, of primary interest. Then one can 
replace the modified Bessel functions in eq. (|37|) by their values at small argument, 
— In(mT) + 0(1), so, to the logarithmic accuracy, one has 

l + 3e-^ 1 



1 — e ^ mT 



The action, eq.(H3|), becomes 



S = Sq\1 — — -7-^ — In ■ 



A 1 — e ^ mT J 

We see that the parameter T enters into the action in the combination A^^ ln(l/mT), 
so in this case it is natural to consider exponentially small T. Let us introduce a new 
parameter, instead of T, 

and assume that 

K ~ 1. 
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These values of T correspond to exponentially large energies: indeed, from the relation 



^ dS 



dT A2T 



we find that 



E 



E 



~ e 



sph 



Analogously, the number of initial particles is 



n 



A fl-e 



-e\2 '^"'sph 



(53) 



(54) 



Since both 6 and k are of order 1, eg. ([5^) implies that the number of initial particles 
is small, 



n 



1 



^sph A 

To calculate the function E, one notices that ET is negligible as compared to S. 
So one has 



E = 2S-ne = 2Sn 



1 



1 + 3e- 



^ 1-e-^ (l-e-^)2 
where 6 and k as functions of E and n are determined by eqs. (|53| ) and 

= arccosh I 1 + %-r — — In ' 



(55) 
Explicitly, 



A2 n 



sph , 



and 



^1 ^ 
K = — m ■ 



According to eqs. (^2]) and (0), the size of the individual instanton increases with 
energy, 

ro oc exp 



, / 1 + 3e-^ ^ 
-A 



1 -e- 

while the separation between the instantons decreases, 

T oc exp(— Aft). 

The validity of the improved dilute instanton gas approximation becomes suddenly 
lost at the energy Ecrit when 



ft ^crit 



1-e- 
3 + 5e- 



(56) 
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It follows from eq. (|55D that at the critical energy, the function F is smaller than 25*0 
only by a factor of order 1, contrary to the case ^ ~ 1/A where the reduction of the 
function F is of order A~^. Thus, at small number of initial particles, n/nsph ~ 
our approximation breaks down when the probability is still exponentially small, with 
the exponent of order of, but smaller than, the instanton one. 

The above formulas simplify considerably in the most interesting case n/usph —>■ 0, 
which occurs in the limit 6 oo. It is clear from eq.(^) that this limit is smooth, 
and the exponent for the probability is 

F{E) = 2So(l-\\n^]. (57) 

Surprisingly enough, eq.(0) coincides with the leading order instanton result, eq. (p7|) . 
Eq.(^) is reliable ok, E < Ecru where Ecru is determined from eq. (|56D to be 

E,r^t ~ e'^/^'Esph. (58) 

This energy is exponentially large in A, but still much smaller than the scale Eq ~ 
e^Esph set by the leading order calculation. The exponent of the probability at ~ 
Ecrit is numerically smaller than the instanton action, 

Fcrit = 2(25*0), 

but the probability is still suppressed exponentially. 

Thus, in the most interesting case of small number of initial particles, the improved 
dilute instanton gas approximation enables one to go beyond the perturbation theory 
about the instanton, but not very far beyond. 



3.5 Multi— instanton amplitudes 

The interest in multi-instanton processes in the context of baryon number violation 
has been raised by the conjecture [PTI 12^ that the multi-instanton contributions 



may become essential at energies where the one-instanton amplitude is still small. 
This conjecture is based on the observation that since the many — > many amplitude 
becomes large at the sphaleron mass, the transition from an initial few-particle state 
that goes through a chain of many many processes (fig. 6) in principle may become 
comparable with the one-instanton one at energies larger than the sphaleron mass. 
This conjecture, that has lead to the premature unitarization hypothesis, is an object 
of controversial discussion. 
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In our model, the function F is reliably calculable at energies above Egph- Fur- 
thermore, as discussed in Sect. 3. 3, the case when n is of order of, but smaller than 
rigph can be treated up to the energy where F becomes parametrically smaller than 
Sq. Thus, our model provides a means to test the relevance of multi-instantons above 
the sphaleron energy. Let us see that the multi-instanton amplitudes are not essen- 
tial when the one-instanton one is small, and, moreover, the p-instanton amplitude 
is exponentially smaller than the {p — l)-instanton one, even at energies larger than 



the sphaleron mass. This will provide a counterexample to the claim of refs.|21, 22 



Consider p-instanton processes at n of order of, but smaller than Ugph- The field 
configuration now consists of 2p instantons located at iT^, j = 1,2, ... ,p and an 
infinite set of fake small instantons at ±Tj ± IT with intensities e"'''^ (the case p = 2 
is shown in fig. 7). For ^ ^ 1, the distances between the instantons are approximately 
equal, so 

^ Ap 

The value of the parameter / is now 

oo p oo 

f = 2j2e-''Ko{lmT)+Y: E e'^^^' MmllT + T, - T,\) 

1=1 j=2l=-oo 
P oo 9 

+ E E e-\'\'Koim\lT + T, + T,\) = p—rM9p/mTp) - ln\ + 0(1), 

j=l /=— oo P 

where Tp and dp may differ from their one-instanton values (it is energy and number 
of initial particles that are fixed). In complete analogy to Sect. 3. 3, one finds for the 
energy and number of initial particles 

while the exponential suppression function is 

Fp = 2pSo (^1 - V^^fi {ep/^Tp)^ ■ (60) 

To calculate the function Fp{E,n) one solves eq.(|59|) with respect to 9p and Tp and 
then substitutes their values into eq. (|60D . From eq.(|59|) it is clear that for given E and 
n, the values of Tp and dp are p times larger than the their values in the one-instanton 
case. So, the p-instanton function Fp{E,n) is p times larger than the one-instanton 
suppression function at the same energy and number of initial particles, 

Fp{E,n)=pF{E,n). 
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We conclude that at those energies when the one-instanton amphtude is exponentially 
suppressed, all multi-instanton contributions are suppressed even stronger. For a 
given number of initial particles, the suppression of multi-instantons disappears at 
the same energy as the suppression of the one-instanton contribution. 

4 Beyond the improved dilute instanton gas ap- 
proximation: periodic instanton 

The improved dilute instanton gas approximation breaks down when the size of the 
core, tq, becomes of order of the distance between the instantons, T, so that the 
instanton cores begin to overlap. At this energy, however, both ro and T are much 
smaller than 1/m, so one may hope to calculate the exponent for the transition 
probability in the leading order in 1/ A beyond the improved dilute gas. Indeed, one 
may expect the exponential term in the lagrangian to play a role only in a small region 
of space-time of order ro or T; in that region the mass term may be neglected, and the 
field obeys the Liouville equation. Outside the non-linearity region, but at distances 
much smaller than 1/m, the field is both free and massless, so that the matching of 
the solutions to the Liouville equation and massive Klein-Gordon equation may still 
be possible. 

It is almost obvious that this procedure should work for periodic instantons at 
all energies up to Egph'- both the improved instanton and the sphaleron have cores 
whose sizes are indeed smaller than 1/m. In this section we construct the periodic 
instantons explicitly in the region of energies where {Egph — E)/Esph ~ 1/A, and the 
improved dilute gas approximation does not work. 

Let us first discuss the region of non-linearity where the field obeys the Liouville 
equation. We do not expect that the solution is 0(2) symmetric (obviously, the 
sphaleron is not), so we have to consider a general solution to the Liouville equation, 




(61) 



where the variable x is introduced in eq. (p!3|) , 



X + ir 



X — IT 



(62) 



z = 



z = 



a 



a 



a is defined in eq.(|T^) and f{z) and g{z) are arbitrary functions. 
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To specify the solution that describes the periodic instanton, we point out first 
that the ordinary instanton is determined by the functions 

fiz) = -, g{z) = -, (63) 

that have poles at the position of the instanton. So, we impose the following require- 
ments on the functions relevant for the periodic instanton: 
i) periodicity, 



ii) existence of poles ai z,z = ilT/2a, Z = 0, ±1, ±2, . . . (recall that in the dilute 
gas approximation, the periodic instanton is a chain of instantons sitting at these 
points); 

iii) regularity of the function x, eg. (pT|) , everywhere in Euclidean space-time; 

iv) reality, 

ensuring that x is real in Euclidean space-time. 

v) reality of x also in Minkowski time (on both Minkowskian parts of the contour 
offig.l). 

These conditions are sufficient to determine / and g. We find 



tanh(/iz) ' tanh(/iz) ' 

where 

27ra 

and B is yet unknown parameter. The field configuration corresponding to this choice 
is 

(64) 



^ 1 



If _ I 

— B cosh nz cosh uz + — sinh uz sinh uz 
/i V B 



The values of T (or yu) and B have to be found by matching this solution to the 
massive free field and requiring that the energy of this solution takes a given value 
E. The contour plots of the configuration (0) for different values of B is presented 
in fig. 8. 

First consider the case B <^1. At small z, z, one has 
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which coincides with eq.(|Tq) ioi c = B / jj. So, B is the diluteness parameter and 
5 -C 1 corresponds to the case of the dilute gas. 

In what follows we concentrate on the case i? ~ 1, when the improved dilute 
gas approximation is not valid. Let us assume that the period T is small enough, 
T <^ 1/m (this assumption will be justified a posteriori). Then the field becomes 
free and massless at T <^ a; <^ 1/m. Indeed, in this region we have fiz ^ 1 and the 
configuration is a time-independent solution to the massless free field equation. 



^ 1 



4/i V B 



\x\ 



We have to match this configuration to the time-independent solution to the massive 
free equation, oc e""^'^'. This matching results in the following relation. 



0(A-^)) 



(65) 



and the field outside the core is 



ve 



-m\x\ 



(66) 



The above formulas explicitly define the periodic instanton at given period T. At 
B = 1, the field (^) is independent of r, and describes the core of the sphaleron. 



A 



In 



/i 



■ cosh fi{z + z) 



The values of the parameter B of order 1 correspond to energies at which the dilute 
instanton gas approximation does not work, i.e., 1 — E/Egph ~ 1/A. The period is 
approximately constant in this energy interval. 



T= = (i + 0(A-^)). 

H Am 



(67) 



Since T does not change much in this region, we make use of the relation OF /BE = 
T/g"^, to obtain that eq.(^) is correct, with the accuracy 1/A^, at all energies including 
the region close to the sphaleron mass. 

To find the value of B, we make use of the following trick. Let us calculate 



explicitly the value of F for the configuration (|6^). We write 

T/2 



F = 2S -FT 



oo 

2 J dt f dx 





^(9.0)^ + 1(9.0)^ + ^0) 



27 



T/2 oo 

2 J dt J dx(^-^{dr(pf + ^{d,(Pf + V{^) 





T/2 oo 

= 2 J dt J dx{dr<p)- 

-oo 

The result of tlie integration for the actual field, eq. 



IS 



A2 



In 



We now compare this result and eq.(^9D and find 



E 



E 



sph 



I- 



This relation determines the value of B at given energy, so that all parameters of the 
solution are now defined. Notice that i? ~ 1 indeed corresponds to {Egph — E)/ Esph ~ 
1/A. Note also that eq.(|57D justifies the assumption that T ^ 1/m, so our calculation 
is rehable at all energies up to Egph- 

The above procedure for obtaining the relevant classical solution is likely to work at 
^ ~ 1/A, i.e., for processes with the number of particles of order ngph- Unfortunately, 
this procedure does not work in the most interesting case ^ ~ 1 (i.e., for small number 
of particles in the initial state): the size of the solution to the Liouville equation very 
rapidly increases with energy and becomes of order 1/m essentially at the same energy 
as one given by eq. ([58|) . The calculations leading to the latter negative conclusion 
are not illuminating and we do not present them here. 



5 Induced vacuum decay 

In this section we briefly discuss a process which is closely related to the shadow 
processes, namely, the decay of the false vacuum induced by initial particles. The 
most probable process of the induced vacuum decay at a given energy and number 
of initial particles is described by the solution to the same boundary value problem, 
eqs.(|D and @, but with the condition that the field in the final state (the part CD 
of the contour of fig.l) describes an expanding bubble with the large value of the 
scalar field inside. In the dilute gas approximation, the configuration is a chain of 
instantons located at (/T, 0) with intensities exp(— (fig. 9). The decay rate is 
again exponential, 

r ~ exp(— F), 
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where F is calculated according to eq.(|D. As is seen from fig. 9, the induced vacuum 
decay can be considered half-instanton process. 

The solution to the boundary value problem can be obtained, in the improved 
dilute instanton gas approximation, in a complete analogy to Sect. 3. The details of 
the calculation are not very instructive, so we present here our results for the rates 
only. 

i) In the case 6 ~ 1/A, the function F{E,n) of the vacuum decay is exactly half 
of the function F for the shadow process at the same energy and the same number 
of particles. The region of validity of the improved dilute gas approximation is also 
the same. 

ii) For parametrically larger ^, 6' ~ 1, the formula analogous to eg. (^51) reads 



So 



e-^ / 9 \ 1^ E 
1 - 2 1 + - In 



1-e-^V 1-e-VA E,ph 
where 6 is again a function of particle number and energy, 

9 = arccosh f 1 H — — — — In ' 



y A2 n EsphJ 

In the limit 9 oo, all instantons, excluding the central one, can be neglected, and 
the function F is equal to its zero-energy value, 5*0. So, we find that at energies where 
the improved dilute gas approximation is valid, the vacuum decay is not enhanced 
by a small (of order f^/A) number of initial particles. The region of validity of the 
dilute gas again extends to a certain critical energy, that depends on the number of 
initial particles and is of order 

Ecrit ~ Q^Esph 

in the limit 9 oo, i.e, when the number of initial particles is small. 



6 (j) model 

In this section we consider another example of the (1+1) dimensional scalar theory 
with the quadratic potential with a cliff. Namely, we take the scalar potential to have 
the following form, 

where is a large number playing the role of A in the potential model. The qualitative 
behavior of the potential is the same as in the exponential model. 
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The sphaleron field outside the core is precisely the same as in the exponential 
model. In particular, the sphaleron energy and the number of particles are given by 
eqs. (pnD and (jl^), respectively. 

The technique developed for the exponential model to describe the instanton and 
improved dilute instanton gas requires a slight modification in the power model. Let 
us find out the instanton in this model. The Euclidean field equation has the form 



c?> = m^0-(iV + 2)^0^+^ (69) 

First let us consider the region r <^ m~^, where the mass term in the field equation 
can be neglected. Though the nonlinear term in eq. (|69|) does not have the exponential 
form, in this region eq.(|69D reduces to the Liouville equation. To see this, let us make 
the following change of variables 



2x 

N 



(70) 



where 0o = 0(0) is a yet undetermined parameter, and consider only such r for which 
x{r) <^ (i.e. 1 — 0(r)/0o ^ !)• Eq.(|69|) obviously reduces to the Liouville equation 
(recall that iV > 1), 

9lx = (71) 



ro 

where 



Nm \ V J 

The 0(2) symmetric solution to eq. (|7TD obeying the condition x(0) = is 

X(r) = In (^1 + . (72) 

We now match the solution, eqs.([70|),([72|), to free massive field at m^^ ^ r ^ ro, 

/ 2 , r2^^ 

and obtain 



r)=0o 1 + ^1 ) =aKo{mr), (73) 



60 = ve'/\ ro ~ m-le-^/^ 
4ei/2y 
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Notice that the size of the instanton core, ro, is indeed exponentially smaller than 
m~^. The instanton action is equal to 



So 



N 

Given the instanton field, it is straightforward to evaluate the leading order in- 
stanton contribution to the two-particle total cross section for the shadow process. 
We find 

... / zL F. \ 

(74) 



Esph / 

SO that the energy scale relevant to the two-particle processes is again exponentially 
large as compared to the sphaleron mass, 

Now let us consider the improved dilute instanton gas. The field outside the 
instanton core is still given by eq.(|32D, while the field inside the core is determined by 



eqs. (|70|) ,(72) where 0o = 0(0) is yet unknown. The analog of the matching condition, 
eq. ([T3|) , is that at mT^ ^ r ^ tq the following relation should hold. 



(r) = 0o(l-^ln + ) =a{Ko{mr)+f), 



where / is defined by the same formula as in the case of the exponential model, 
eq.(|38|). One obtains 

(po = v exp Q - , ro ~ m"^ exp (^-^ + . 

It is worth noting that the parameter a depends on /, i.e., on the parameters T and 
6, in contrast to the exponential model, where a is constant. So, not only inside 
the core, but also outside the core the field is sensitive to the distance between the 
instantons. 

In the case 9 ~ i.e., at n/risph ~ 1, one obtains from eqs.(^5D, (pif ), (^) and 
( |75D the following expressions for the energy and the number of initial particles, 

E = x'exp{l-xf\if3))f2{f3)Esph, 
n = x^ exp(l - xfi{P))f3{P)nsph, 
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where 

Svr 

^ ~ NmT 



/3 



mT 

and /i, /2 and /s are the same functions as defined in Sect. 3. The calculations 
analogous to those of the exponential model give for the action 

^=^exp(l-x/i(/3)), 
so the function F = 2S — FT — nO has the following form, 

:i-x/i(/3))exp(l-x/i (/?)). 



Though these formulas are slightly more complicated than eqs.(^6D and (0), the 
qualitative behavior of the function F{F, n) is the same as in the case of the expo- 
nential model. In other words, the behavior shown in fig.5 is qualitatively correct in 
the power model as well. The improved instanton gas approximation is valid until 
the probability becomes parametrically less suppressed as compared to the instanton, 

iV' 

in a complete analogy to the exponential model. 

In the regime 5 ~ 1, i.e., at n/usph ~ 1/^5 interesting energies are exponen- 
tially large, so that 

— m — — = K ~ 1. 

N Esph 



The formula analogous to eq.(^5|) is 



16ge-^ ■ 
(l-e-^)2'^ 



F = 2Soexpi-A^-—^K 

where 6 is to be determined from a relation analogous to eq.(^ 

QATT^eUsph e-^ ( 1 + 3e-^ ^ 



The size of the core is 



1 / N ^ l+3e-^ ^ 

ro ~ — exp -— + A^- —k 

m V 4 1 — e ; 
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In the limit 9 oo, i.e., n/nsph — ^ 0, the behavior of the function F is as follows, 

Note that in the 0^ model, the function F{E/ Esph) contains higher order terms in 
N~^ln{E/Esph), in contrast to the exponential model, where all higher terms vanish 
and the result coincides with the one-instanton expression. 

Let us explain the latter result. Recall that in the exponential model only the 
field inside the instanton core is modified by the presence of other instantons, while 
the linear tail of an instanton in the chain is the same as that of a single instanton 
(the parameter a is constant independent of T and 9). If the cores do not overlap, 
the field configuration outside the cores is the same as in the ordinary dilute instan- 



ton gas approximation. The dependence of E and n on T and 9, eqs.([3^) and (pSf), 
thus coincides with the formulas of the ordinary dilute instanton gas approximation. 
Since the function F can be recovered from eq.(P), it is not surprising that the im- 
proved dilute instanton gas implies the same result for F as the lowest order of the 
perturbation theory around a single instanton. In the 0^ model, on the contrary, the 
intensity of the linear tail depends on T and 9 (see eq.([75|)), so the correction to the 
one-instanton formula for the function F is large. 

In the limit 9 ^ oo, the improved dilute gas approximation breaks downs at 
Ecrit ~ e'^^^Egph, where F = 2So/y/e. As in the exponential model, the function F is 
reduced only by a factor of order one ( 2/3 in the exponential model and e~^/^ in the 
power model). 



7 Conclusions 

In this paper we have found the solution to the classical boundary value problem 
for certain (1 + 1) dimensional field theories with the potential having the form of a 
quadratic potential with a cliff. The presence of an additional large parameter in 
the model allowed us to deal with energies comparable to or much larger than the 
sphaleron mass. The solution, which describes the shadow processes, has been found 
reliably within the improved dilute instanton gas approximation at large energies 
where the ordinary dilute instanton gas approximation does not work. We have found 
that for the number of initial particles, n, of order ngph, the function F{E/Esph,n) 
monotonously decreases from 25*0 at _E = and becomes of order Sq/X (or Sq/N) 
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at some Ecru, where the improved dilute gas approximation fails. For n of order 
nsph/^ (or Usph/N in the (f)^ model), including the limit n ^ which is smooth, the 
improved dilute gas approximation breaks down at energies where the function F is 
still of order 5*0, and the instanton-like transitions are still suppressed. So, although in 
our models we can go beyond the ordinary perturbative theory around the instanton, 
the possibility of large probability of the instanton-like transitions induced by two 
energetic particles remains an open question even in these models. 

The authors are indebted to A.N.Kuznetsov and P.G.Tinyakov for numerous help- 
ful discussions at all stages of this work. We are indebted to T.Banks and V.Petrov 
for the discussion of the results. V.R. thanks Aspen Center for Physics and Rutgers 
University for hospitality. The work of D.T.S. is supported in part by the Rus- 
sian Foundation for Fundamental Research (project 93-02-3812) and by the Wein- 
gart Foundation through a cooperative agreement with the Department of Physics at 
UCLA. 
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Figure captions: 

1. The contour in the complex time plane, on which the boundary value problem 
is formulated. 

2. The potential V{(j)). 

3. The schematic plot of the solution to the boundary value problem in the 
improved dilute instanton gas approximation. / stand for the instanton. 

4. The function F{E) for the periodic instanton. 

5. The function F{E,n) for u = n/usph ~ 1- In the regions where F ~ Sq/X, the 
curves are unreliable. 

6. Multi-instanton contributions to 2 ^many amplitudes. 

7. The two-instanton solution to the boundary value problem. 

8. The contour plots of the periodic instanton solution, eg. (p^) , for differnent 
values of B. Larger B correspond to higher energies. As B increases from {E = 0) 
to 1 {E = Egph), the configuration evolves from a chain of instantons to the sphaleron. 

9. The configuration decribing the induced vacuum decay in the improved dilute 
instanton gas approximation. 
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